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Abstract 

We study the dynamics of a self-gravitating cooling filamentary cloud using a simplified model. We 
concentrate on the radial distribution and restrict ourselves to quasi-hydrostatic, cylindrically symmetric 
cooling flows. For a power-law dependence of cooling function on the temperature, self-similar solutions 
which describe quasi-hydrostatic cooling fiows are derived. We consider obtically thin filaments with a con- 
stant mass per unit length and the solutions are parameterized by their line masses. There is no polytropic 
relation between the density and the pressure. The filament experiences radiative condensation, irrespective 
of the 7, the gas specific heat ratio. So, the filament becomes denser due to the quasi-hydrostatic fiows and 
the density at the center (pc) increases in proportion to {to — t)~^ , where t denotes the time. The term, 
to, denotes an epoch at which the central density increases infinitely. We also found that the radius of the 
filament (rc) decreases in proportion to (to — t)'^ ■ 

Subject headings: hydrodynamic - instabilities - ISM: clouds - stars: formation 

1 INTRODUCTION 

Interstellar molecular clouds are very complex structures in which many physical processes come into play. 
In these clouds, rich and complex structures have been seen (see, e.g., Myers 1991). Since most star forma- 
tion essentially occurs in molecular clouds, understanding the origin and evolution of their structures is very 
important. 

Filamentary structures associated with clumps and cores are very common (e.g., Schneider & Elmegreen 
1979; Houlahan & Scalo 1992; Wiseman & Adams 1994; Alves et al. 1998; Harjunpaa et al. 1999). It is now 
widely accepted that this hierarchical structure has been formed by gravitational instability. In this scenario for 
the formation of the hierarchical structure, four distinct phases have been studied: (a) gravitational instability 
along the filamentary cloud leads to the formation of clumps (e.g., Nagasawa 1987; Nakamura, Hanawa, & 
Nakano 1993; Gehman, Adams, & Watkins 1996; Nagai, Inutsuka, & Miyama 1998; Fiege & Pudritz 2000); 
(b) these clumps become geometrically thin disks, perpendicular to the filament axis (Nakamura, Hanawa, & 
Nakano 1995; Nakamura et al. 1999); (c) The disks collapse toward their centers and, because of gravitational 
instability, each of them fragments into cloud cores (Nakamura & Hanawa 1997; Saigo & Hanawa 1998); and 
finally (d) these cores are scattered by gravitational interaction. 

It must be noted that in this scenario, the overall collapse of a filamentary cloud is slower than the growth 
of local density fluctuations. In fact, it has been assumed that the parent filamentary cloud is well below the 
critical mass per unit length required for radial collapse to a spindle. There is a possibility that the filamentary 
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cloud collapses, however, and several authors have studied the collapse of filaments (Bonnell & Bastien 1991; 
Inutsuka & Miyama 1993; Kawachi & Hanawa 1998). 

Most of the previous works have studied the collapse of spherically symmetric clouds (e.g., Larson 1972; 
Terebcy, Shu, & Casscn 1984; Foster & Chevalier 1993). Calculations showed that the temperature of a 
collapsing spherical molecular cloud remains roughly constant in the early phases of star formation (see, e.g., 
Larson 1969). Investigating the collapse of a spherical cloud using three dimensional numerical simulations is 
very difficult. But there exists asymptotic self-similar solutions for the collapse of a spherical cloud (Larson 1969; 
Penston 1969; Shu 1977). Also, self-similar solutions for the collapse of a filamentary cloud were investigated, 
and different sets have been found (Inutsuka & Miyama 1992; Kawachi & Hanawa 1998; Semelin, Sanchez, & 
dc Vega 1999). Kawachi & Hanawa (1998) investigated the gravitational collapse of a filamentary cloud using 
zooming coordinates (Bouquet et al. 1985). They used a polytropic equation of state to indicate the effects of 
deviations from isothermality in the collapse. 

In fact, all of these authors followed the usual (Jeans) scenario for star formation. But as we shall show 
in this paper, there exists another (non- Jeans) scenario for star formation in filamentary clouds. Because of 
radiative cooling, a filamentary cloud which is in equilibrium, may evolve over a longer period of time. Thus, 
radiative condensation flows can be formed in a cooling filamentary cloud. So, although the cloud is in a 
quasi-hydrostatic state, we may have density collapse, due to the presence of a radiative condensation flow. 
The same problem has been investigated for the dynamics of cooling spherical clouds (Meerson, Megged, & 
Tajima 1996; hereafter MMT). In the present work, we have studied radiative cooling flows in a filamentary 
cloud, using a family of similarity solutions and instead of a polytropic equation of state, the energy equation 
is used. We found solutions that did not indicate any polytropic relation between pressure and density. Since 
we have considered the constraint of a constant mass per unit length, the solutions were parameterized by their 
line masses. We show that for the power-law radiative cooling function, the dynamics of quasi-hydrostatic, 
cylindrically symmetric flows in a filamentary cloud is independent of the specific heat ratio of the gas. 

The plan of this paper is as follows. We present the general mathematical framework and equations in § 2. 
By introducing similarity solutions, the dynamics of a radiative cooling filamentary cloud has been studied in § 
3. We concentrate on quasi-hydrostatic flows. In § 4, after discussing the behaviour of central density and the 
radius of the cloud, we conclude with the implications of our model. 



2 GENERAL FORMULATION 

We present here the basic equations used to describe cylindrically symmetric radiative cooling flows. In this 
work, we neglect the effects of viscosity and magnetic fields. In cylindrical coordinates {r,(fi,z), which the 
filamentary cloud is assumed to be long in the 2— direction, we have: 

dp 1 d , . ^ 

,dv dv .^ dp d'^ , , 

19 a* 

-j-{v—)=^'kGp, (3) 
r or or 

where p,v,p,'^ denote the gas density, radial velocity, pressure, and gravitational potential, respectively. We 
write the energy equation and the equation of state as follows: 

P=^l/T. (5) 

where A{p,T) is the radiative loss function and Rg and iig are the gas constant and the effective molar mass, 
respectively. Neglecting any external heating, we can fit the radiative loss function by a power law. 



A(p,T)=A.p2T^ (6) 
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where v and Ay are constants. Depending on selected intervals of temperature, these parameters can be 
determined (Spitzer 1978). 

Our scope of interest is quasi-hydrostatic flows in filamentary clouds. In fact, this situation occurs when 
cooling is slow on the dynamic time scale. For quasi-hydrostatic flows, the equation (2) becomes 

The same relation has been used by MMT. 

To simplify the problem, we introduce the dimensionless variables according to 

p —>■ pp,p pp, * -» **,T -» fT,v —>-vv,r—>- rr,t -» it, (8) 

where 



po,p = po,f = To,^ = P^,f = {-^)i,i=-, 2rru ^^ = T (9) 

Po 47rGp^' (7-l)A^p^To^ t 



Under the transformation equations (8), the equations (1) and (7) do not change and the rest of the equations 
are cast into these forms: 

3 SELF-SIMILAR SOLUTIONS 

It is useful to study the similarity solutions of equations (1), (7), (10) and (11). We introduce the similarity 
variable as y = where ro{t) = {to — t)^ and t <to- Using this transformation, we have the following forms 

of p, p, V and 

p{r,t) = {to-t)-^R{y), (12) 
p{r,t) = {to-t)-^P{y), (13) 
v{r,t) = {to-tr"-V{y), (14) 
*(r,i)=S(y), (15) 
where R{y)., P{y), V{y) and S{y) satisfy these equations: 

If - (") 

dP dS 

y dy a.y 
11 dP P d 

P+{^+V)—+^--{yV)+R^-''P'' = Q. (19) 
2 dy y dy 

This form of similarity solutions deals with blow-up, or at least rapid growth, in a finite time to- Since we 
arc interested in cooling fiows of an isolated filamentary cloud, the constraint of a constant mass per unit length 

can be written as 



prdr = m = const, (20) 
where m = ^^^M and M is mass per unit length. 
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Now, wc have a set of ordinary differential equations. It is interesting that equations (16), (17), (18) and 
(19) are integrable. Using equation (16) we obtain 

yRV+'^y^R = c, (21) 

where c is a constant. The velocity of the flows on the axis of a filamentary cloud is assumed to be zero. So we 
obtain 

V{y) = -|. (22) 

It shows that the velocity of the cylindrically symmetric cooling flow is independent of u, the exponent of the 
power-law cooling function. MMT showed that for spherically symmetric cooling flows, we may have gas inflow 
or outflow, depending on the value of v . Substituting equation (22) into (19), we find 

P(y) = (7-l)"[i?(y)f, (23) 

where a = and (3 = This relation shows that v may have any positive or negative values except v = 1 

or 2. In addition there is no polytropic relationship between the pressure and the density. 
Finally, after substituting equation (23) into equations (17) and (18) we obtain 

where 

e = [i?(j/)]i,^ = (n+ l)-5(7 - l)^t/,n = 1 - i/. (25) 

Equation (24) represents the form of the Lane-Emden equation appropriate for cylindrical polytropes. This 
equation on the study of the equilibrium of gaseous filaments has been studied (Ostriker 1964). The boundary 
conditions are 

rift 

^(^ = 0) = 1,(^)^=0 = 0. (26) 

It can be easily shown that solutions of the cylindrical Lane-Emden equation satisfying boundary conditions 
(26) miist decrease monotonically from the origin and eventually approach zero. For n = (i^ = 1) or n = 1 
(i/ = 0), this equation becomes linear and the solutions can be readily obtained. If 9n{C) reperesnts the required 
solution, we have 

Oo{^) = i-\eMO = MO- (27) 

We note that in our problem, ^o(0 '^o* a solution. For other values of n, the solutions are not in closed forms. 
In general, however, Ostriker (1964) obtained a series solution: 

- 1 - ^ + " .4 _ ^(3^ -2) ,6 , »il8/>^-29» + 12) ^3 _ 

(1!)222^ (21)22^^ (3!)226 ^ ^ (4!)228 ^ ^ ' 

We note that the solutions must satisfy the constraint of conservation of mass. In fact, the solutions of 
the cylindrical Lane-Emden equation are subject to the same homology transformations as are the solutions of 
the spherical Lane-Emden equation: if 6n{C) is a solution, then A"^0„(^^) will also be a solution (n ^ 1). 
Using this transformation, it is possible after finding the solution, then rescale it. Therefore, equation (20) is 
as follows: 

t\9n{A^)nd^ = ^(7-1)-, (29) 

where is the first zero of 9n{Cl- AH of the solutions 6'„(^) have a first zero; i.e., the radii and masses of the 
filamentary clouds are finite. For mathematical proof of this point see Ostriker (1964). If n = 1 {v = 0), then 
A6\{^) will be a solution and the equation (29) is replaced by 

Aj%,{0^d^ = ^{^-1), (30) 
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where ^i(^) = Jo(0 , so 

Equation (12) shows that for quasi-hydrostatic cylindrically symmetric coohng flows, the density at the 
center (pc) increases in proportion to {to — t)^^. This is slower than the growth of density at the center, 
during gravitational collapse, which is in proportion to {to — t)~^ (Kawachi & Hanawa 1998). We see that in 
a filamentary cloud Pc is independent of i^, while in a spherical cloud the central density Pc depends on v and 
varies in proportion to {to — t)^2+77 (MMT). 

Since the radius of the cloud (rc) is defined by the first zero of 0n{S,), we can see that Vc decreases in 
proportion to {to — t)^, which is independent of v. MMT showed that in spherical clouds, Tc is proportional to 
{to — t) which depends on v. 

Another important difference between cylindrical and spherical cooling flows exists. In fact, MMT showed 
that the dynamics of quasi-hydrostatic spherically symmetric cooling flows depends on 7, the gas specific ratio. 
They found that 7 = | is a critical value, i.e., if 7 is greater (smaller) than |, then the cloud condenses 
(expands). But our study showed that in quasi- hydrostatic cylindrically symmetric cooling flows, there is no 
critical value. In fact, for all values of 7 > 1, filamentary clouds undergo radiative condensation. 



4 DISCUSSION AND CONCLUSIONS 

In this paper, we have studied quasi-hydrostatic radiative cooling flows of self-gravitating filamentary clouds. 
A non- Jeans scenario for star formation in these clouds has been investigated. This scenario in spherical clouds 
was presented, as far as we know, for the first time by MMT. 

When the line mass of a filamentary cloud exceeds the critical line mass, the entire filament collapses towards 
the axis without fragmentation. But a filamentary cloud which is in equilibrium is unstable along its axis and 
it will fragment. In this paper we suggested another possibility to this view point. A filamentary cloud which 
is in equilibrium may evolve over a longer time-scale, i.e., due to the cooling of the cloud, quasi- hydrostatic 
flows can be formed. Considering a simple power-law cooling function, we showed that the filament experiences 
radiative condensation. These cooling flows can be parameterized by the line mass of the filament. 

It is remarkable that the decrement of the radius and the growth of density at the center of the cloud are 
independent of v: i.e., we have rc oc {to — t)^ and pc oc {to — t)~^. We must note that these results strongly 
depend on the form of the cooling function. But it seems that we can use piece- wise solutions. That is, consider 
that for some part of the contraction, A„ and v have one set of values, and then join this on to another phase 
of contraction, when Ai, and v have another set of values. What will change is to, but the equations should join 
with initial and boundary conditions in this way, it seems. If this is true, then we can, in principle, consider 
a much wider range of cooling functions with the same analytical-type solution. By finding suitable values of 
Ai, and for the range of densities and temperatures where there is a good approximation, we can find the 
analytical solution advancing forward from the final state of the previous partial solution. Although we actually 
do not do this in this paper, using this method we can broaden the applicability of our analysis. Our results 
woiild also apply to a turbulent filament, for some range of parameters. In fact, our cooling function in equation 
(6) is a general function and it might include, in the same qualitative sense, turbulent energy dissipation as 
well. In other words, one could represent the dissipation rate of turbulent energy as a power law dependent 
on density and velocity dispersion, and identify the velocity dispersion with the temperature in our equations. 
There are some nice studies of turbulent dissipation (see, e.g., Mac Low 1999). We suggest that, at least, for 
some range of parameters, our cooling expression might be used as an approximation to the turbulent cooling 
rate. Then our results would also apply to a turbulent filament for that parameter range. 

Since the behaviours of rc and pc obtained based on the similarity solutions, it must be confirmed by a 
numerical calculation of solutions. We can see that for other forms of cooling function, finding self-similar 
solution is much more difficult. Thus, numerical simulations will reveal the behaviour of cooling flows in 
filaments. Accordingly, initial conditions which will lead to the formation of such flows can be determined. We 
expect that the similarity solutions for describing quasi-hydrostatic cooling flows in fllamentary clouds will serve 
as guidelines for understanding fllaments and doing numerical simulations. We must note that our solutions 
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apply only to the early stages of filament contraction, and after the density in the center becomes high enough, 
perhaps the Jeans-collapse solutions would take over. Even if the Jeans-collapse is eventually important, we 
have still gained valuable insight into the early stages of filament formation. Since cosmological filaments, 
diffuse HI filaments, and molecular filaments, all of which show some cooling properties, our solutions might 
have applications to the formation of these structures. 

Each molecular filament typically may contain several distinct clumps which are spaced along that filament. 
We can consider the filament and the embedded clumps as a self-gravitating system which is in equilibrium 
(Curry 2000). Theoretical arguments suggest that the mass of each clump is much larger than a typical stellar 
mass (see, e.g., Fiege & Pudritz 2000). In our scenario, on a longer time-scale the mass and the shape of these 
clumps have been affected by quasi- hydrostatic cooling flows. It will be interesting to investigate the effects of 
cooling flows on the clumps. 
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